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£S) , Abstract 

^H ■ 

£> ' The production rate of prompt ij/'s at large transverse momentum at the Tevatron 

\Q . is larger than theoretical expectations by about a factor of 30. As a solution to this 

m ■ , 

puzzle, we suggest that the dominant ij) production mechanism is the fragmentation of 

t^J- , a gluon into a cc pair in a pointlike color-octet S-wave state, which subsequently evolves 

ON , 

^5 | nonperturbatively into a ip plus light hadrons. The contribution to the fragmentation 

function from this process is enhanced by a short-distance factor of 1/a^ relative to the 

conventional color-singlet contribution. This may compensate for the suppression by 



v 4 , where v is the relative momentum of the charm quark in the tp' . If this is indeed the 



^/ dominant production mechanism at large pr, then the prompt ip s that are observed 

at the Tevatron should almost always be associated with a jet of light hadrons. 



The production of heavy quarkonium in high energy processes provides severe tests 
of our understanding of QCD in a domain that should be accessible to theoretical analysis. 
At high energies, QCD mechanisms for the direct production of charmonium are embedded 
in a large background from the production of b quarks, followed by the decay of the resulting 
B hadrons into charmonium states. Vertex detectors can be used to separate the charmo- 
nium states coming from b quarks from those that are produced by QCD interactions, thus 
allowing the direct study of the prompt QCD production mechanisms. In the 1991-92 run 
of the Tevatron, the CDF detector group used their new vertex detector to study prompt 
charmonium production |l], |2|. The results have proved to be an embarrassment to theory, 
revealing orders-of-magnitude discrepancies with theoretical predictions. 

Until recently, the conventional wisdom was that the dominant mechanisms for prompt 
charmonium production at large transverse momentum px come from the processes such as 
gluon-gluon fusion which are leading order in the QCD coupling constant a s 0. The result- 
ing predictions for prompt J/ip production at large px are smaller than the CDF data by 
more than an order of magnitude. However, as pointed out by Braaten and Yuan in 1990 M, 
the dominant production mechanism at sufficiently large pt is the production of a parton 
with large transverse momentum, followed by the fragmentation of the parton into char- 
monium states. Including this fragmentation mechanism brings the theoretical predictions 
for prompt J/ip production to within a factor of 3 of the data || |6| , close enough that the 
remaining discrepancy can be attributed to theoretical uncertainties. Unfortunately, even 
after including the fragmentation contribution, the prediction for the %p' production rate 
remains a factor of 30 below the data. This remarkably large discrepancy between theory 
and experiment suggests that an entirely new production mechanism must be dominating 
the production of ip' at large px- 

Previous proposals for solving the ip' surplus problem have focused on the production 
of higher charmonium states which decay into ip' . Among the possible states are higher P- 
wave states, higher D-wave states, and states whose dominant Fock component is \ccg) @. 
For each of these possibilities, there is reason to question whether the production rate can be 
large enough and whether the branching fraction into ip' can be large enough to account for 
the observed rate for prompt ip' production. In this Letter, we propose a completely different 



solution to the if)' surplus problem. We suggest that the dominant production mechanism for 
ip' at large p? is gluon fragmentation into a cc pair in a pointlike color-octet S-wave state, 
followed by the nonperturbative evolution of this state into a ip' plus light hadrons. We 
calculate the corresponding contribution to the fragmentation function for g — > ip' in terms 
of a well-defined nonperturbative matrix element, and we determine the value of this matrix 
element that would be required to explain the production rate of ip' at the Tevatron. We 
then discuss possible tests of this solution to the ip' surplus problem. 

Our proposal is based on a general factorization analysis of the annihilation and 
production of heavy quarkonium by Bodwin, Braaten, and Lepage ||. The factorization 
formalism of Ref. || allows the fragmentation functions for charmonium production to be 
factored into short- distance coefficients, which describe the production rate of a cc pair within 
a region of size l/m c , and long-distance factors that contain all the nonperturbative dynamics 
of the formation of a bound state containing the cc pair. The fragmentation function for a 
gluon to split into a charmonium state X with longitudinal momentum fraction z can be 
written 

D^ x {z,n) = Y,dn(z,»)(0\OX\0), (1) 

n 

where O* are local 4-fermon operators that were defined in Ref. || in terms of the fields 
of nonrelativistic QCD ||. The short- distance coefficients d n (z,fi) are independent of the 
quarkonium state X. For a fragmentation scale fi of order m c , they can be computed using 
perturbation theory in a s (m c ). The dependence on the quarkonium state X appears in the 
long-distance factors (O^). The relative importance of the various matrix elements for a 
given state X can be estimated by how they scale with v, the typical relative velocity of the 
charm quark in charmonium M. 

In previous calculations of fragmentation functions for charmonium production, only 
the matrix elements that dominate in the nonrelativistic limit v — > have been considered. 
In the case of ip', the only such matrix element in the notation of Ref. || is (Of ( 3 Si)). In 

this vacuum matrix element, the local operator Of ( 3 Si) creates a cc pair in a color-singlet 
3 Si state, projects onto the subspace of states which in the asymptotic future contain a if)', 
and then annihilates the cc at the creation point. Up to corrections of relative order v 4 , this 



matrix element can be related to the wavefunction R^i /v47r of the if)' at the origin: 



T0V3, 



(0|On d 5i)|0> « -|^|- (2) 

In the case of S-waves, the factorization approach reduces in the nonrelativistic limit to the 
color-singlet model fllO |, where the if)' is treated as a cc pair in a color-singlet 3 S± state with 



vanishing relative momentum. 

In the case of P-waves like the Xcj state, there are 2 matrix elements that survive 
in the nonrelativistic limit. In addition to (Oi cJ ( 3 Pj)) , which is related to the derivative 
of the wavefunction at the origin, there is also the color-octet matrix element (Og CJ ( 3 Si)) . 
The local operator Og CJ ( 3 Si) creates a cc pair in a color-octet 3 Si state, projects onto the 
subspace of states which in the asymptotic future contain a Xcj, and then annihilates the 
cc at the creation point. In the color-singlet model, only the matrix element ((9* cJ ( 3 Pj)) is 
taken into account. This model is perturbatively inconsistent because of infrared divergences 
in the radiative corrections to P-wave production rates. In the factorization formalism of 
Ref. ||, this problem is solved by factoring the infrared divergences into the color-octet 
matrix element (£$ cJ ( 3 Si)) [[II]. 



The restriction to the leading matrix elements in the v — > limit is a great simplifi- 
cation, but it is not necessarily correct. The typical value of v 2 in the case of charmonium 
is only about 1/3, so factors of v 2 do not provide a very large suppression. In the case 
of ([!]), dropping contributions that are suppressed by powers of v 2 might be reasonable if 
the short-distance factors d n (z,[i) were all comparable in size. However these factors are of 
different order in a s for different matrix elements. Short-distance factors of a s (m c )/7r prob- 
ably provide much more effective suppression than factors of v 2 . Thus, in estimating the 
importance of various terms in the fragmentation function, one should take into account not 
only the scaling of the matrix element with v 2 , but also the order in a s of the short- distance 
coefficient. 

In the case of gluon fragmentation into if)', the short-distance factor associated with 
the matrix element (Of ( 3 5i)) is of order a 3 . It describes the process g* — > ccgg, with the cc 
pair in a color-singlet 3 S\ state. By the velocity-scaling rules of Ref. H, the matrix element 



(Of ( 3 Si)) is of order m 3 t> 3 , so the contribution to the fragmentation function is of order 
a 3 v 3 . While all other matrix elements are suppressed by powers of v 2 , there are some with 
fewer short-distance suppression factors of a s . In particular, there are some matrix elements 
for which the leading-order short-distance process is g* — > cc and the rate is of order a s . 

Of these matrix elements, the leading one in the nonrelativistic limit is (Of ( 3 Si)), which 
scales like m 3 t> 7 . The suppression factor of t> 4 relative to (Of ( 3 Si)) reflects the fact that 
the cc pair in the ip' can make a double El transition with amplitude of order v 2 from the 
dominant \cc) Fock state to a \ccgg) state in which the cc pair is in a color-octet 3 5*i state. 

The contribution to the fragmentation function for g — > ip' from the short- distance 
production of a cc pair in a color-octet 3 S\ state was calculated to leading order in a s in 
Ref. ||12| . The calculation was carried out for the specific case of c/ — 5. ^ cJ , but since the 
short-distance coefficients in ([I]) are independent of the quarkonium state, the result applies 
equally well to ip' . The fragmentation function is 

JW(*,A«) = ^^^ l ~ z ) <0|<?f( 3 Si)|0>. (3) 

In Ref. [0, the corresponding fragmentation function for g — > Xcj was expressed in terms 
of a quantity H' s defined by 

H > - wh< (m " {3sm - (4) 

The form (|3|) is preferred, since it corresponds to the factorization of the fragmentation func- 
tion into a short-distance factor proportional to a s jrn? c and a long-distance matrix element. 

In the fragmentation function for g — + Xcj, the matrix element (Og cJ ( 3 Si)) can be 
determined phenomenologically from data on charmonium production in 5-meson decays. 
From eq. (24) of Ref. []TTJ , we have 

T(B^ Xc2 + X) _ liO j0\Or 2 ( 3 Si)\0) () 

T(B -»• ev e + X) ' m 2 c m b ' l ' 

The branching fraction for B — ► \c2 + X has recently been measured by the CLEO collabo- 
ration to be (0.25 ±0.11)% [[H|. Assuming m c ~ 1.5 GeV and m& ~ 4.5 GeV, the expression 



(0) gives 

(0\Og*( 3 Si)\0) = (0.016 ± 0.007) GeV 3 , (6) 

which corresponds to H' & rj (1.4 ±0.6) MeV. This value is significantly smaller than the esti- 
mate H g « 3 MeV obtained in Ref . [[lip and used in Refs. and . Thus the contribution 
to J/tp production from gluon fragmentation into Xcj that was calculated in Refs. || and 
[[J should be decreased by a factor of 2. The theoretical prediction for J/ip production then 
falls below the data by about a factor of 5. 

In the case of ip', we can determine the matrix element (Of ( 3 Si)) by fitting the CDF 
data on the production rate of prompt if}' at large pr- The production rate is given by 
convoluting hard scattering cross sections for the production of gluons with large transverse 
momentum with parton distributions for the colliding proton and antiproton and with the 
fragmentation function D g _,^i(z, n). The fragmentation function given by (§) is evolved up 
to the scale \x = pr/ z of the transverse momentum of the fragmenting gluon by using the 
Altarelli-Parisi equations. The production rate measured by CDF || in the 1991-92 run is 
shown in Figure 1, along with three theoretical curves. The dotted line is the production rate 
from leading-order production mechanisms, such as gg — > ip'g, while the dashed line is the 
color-singlet fragmentation contribution calculated in Ref. ||, which is dominated by charm 
quark fragmentation into if)'. The solid curve is the color-octet fragmentation contribution, 

with the matrix element (Of ( 3 Si)) adjusted to make the normalization agree with the data. 
The resulting value of the matrix element is 

(0|C»f( 3 5i)|0) = 0.0042 GeV 3 . (7) 

Note that the shape of the color-octet fragmentation contribution fits the data perfectly. 

There is a very simple consistency check that the color-octet fragmentation mechanism 
must satisfy in order to provide a viable explanation for the if}' surplus. The value of the 

matrix element (Of ( 3 Si)) must be small enough relative to the matrix element (Of ( 3 5'i)) 
to be consistent with suppression by a factor of v 4 . The color-singlet matrix element can be 
determined from the leptonic decay rate of the ip' to be 0.11 GeV . The color-octet matrix 
element (|7|) is smaller by a factor of 25, consistent with suppression by v A . 



There are several experimental consistency checks of our proposed explanation of the 
ip' surplus. If color-octet fragmentation is indeed the dominant production mechanism, then 
most of the if> n s produced at the Tevatron must be accompanied by light hadrons. The 
leading order mechanisms, such as gg — > ijj'g, produce isolated if) n s with no accompanying 
hadrons. Color-singlet fragmentation produces ^"s that are associated with a charm quark 
jet, but with transverse momentum of order m c relative to the jet. Color-octet fragmenta- 
tion also produces ^"s that are associated with a jet of light hadrons, but the transverse 
momentum relative to the jet is of order m c v or smaller. Thus measurements of the hadronic 
energy distribution as a function of the distance from the if)' momentum in azimuthal angle 
and in rapidity could provide evidence for this new production mechanism. It has recently 



been pointed out by Cho and Wise |14| that measurements of the spin alignment of the 
if>"s, which is reflected in the angular distribution of their leptonic decays, might provide a 
clear signature for the color-octet fragmentation mechanism. This mechanism produces if) n s 
that are 100% transversely polarized at leading order in a s , while other mechanisms tend to 
produce unpolarized ip n s. 

If color-octet fragmentation is so important in the case of ip' production, one might 
ask why it is not equally important for J/ip production at large pt- In fact, it might very well 
be. The main difference is that the J/ip signal is fed by gluon fragmentation into Xci and Xc2 
followed by their radiative decays into J/ip, and this overwhelms the rate for color-singlet 
fragmentation directly into J /ip. However there is room for a substantial contribution from 
color-octet fragmentation into J/ip, especially with the new determination @ of the matrix 
element (Oq c2 ( 3 Si)), which leaves the predictions of Refs. || and || a factor of 5 below the 

CDF data. If the matrix element (O s ( 3 Si)) for J/ip has a value comparable to that for 
ip' given in (0), then the color-octet fragmentation mechanism for J/ip production will be 
comparable in importance to gluon fragmentation into %& an d Xci- 

If a color-octet contribution dominates the fragmentation function for g — > if)', one 
might worry that annihilation of a cc pair in a color-octet state should also dominate the 
decay rate of the ip' into light hadrons. However the relative importance of the color- 
singlet and color-octet contributions is different in the case of annihilation. The color-singlet 
contribution to the decay rate is of order a 3 m c v 3 , while annihilation from a color-octet S-wave 



state contributes at order a 2 s m c v 7 . Thus the short- distance enhancement of the color-octet 
contribution is only one power of l/a s , rather than 2 powers as in the case of fragmentation. 
Furthermore, there is no apparent rigorous relation between the color-octet annihilation 
matrix elements and the color-octet production matrix elements beyond the fact that they 
have the same scaling in v M. 

Our proposal for the solution of the if/ surplus problem has predictive power, because 
the matrix element (Of ( 3 Si)) determined from ip' production at large px can be used to 
predict ip' production in other high energy processes. Unfortunately, the matrix elements 

(Of ( 1 5'o)) and (Of ( 3 Pj)) are of the same order in v 2 as (Of ( 3 Si)), and they must also 
be taken into account. Different production processes probe different linear combinations of 
these 3 color-octet matrix elements at leading order in a s . Gluon fragmentation is rather 

unique in that (Of ( 3 Si)) is the only color-octet matrix element that contributes at leading 
order in a s and in v 2 . In photoproduction of ijj', on the other hand, only the matrix elements 
(Of ( 1 5'o)) and (Of ( 3 Pj)) contribute at leading order in a s and in v 2 . While suppressed 
relative to the color-singlet contribution by a factor of v 4 , they are enhanced by a short- 
distance factor of l/a s . Moreover, at leading order in a s , the color-octet contribution is 
concentrated near z = 1, where z is the ratio of the energy of the ip f to that of the photon. The 
cross section is known experimentally to peak strongly in the z — > 1 region, and charmonium 
states produced in this region are labelled "elastic". Thus a linear combination of the 

matrix elements (Of ( 1 5'o)) an d (Of ( 3 Pj)) can be determined from the elastic peak in ip' 
photoproduction 111 



In most of the work on the production of charmonium over the last decade, the 
dominant production mechanisms have been assumed to be ones in which the cc pair is 
produced at short distances in a color-singlet state. However, in the case of P-wave states, it is 
essential to include the color-octet mechanism to obtain consistent perturbative predictions. 
In the case of S- waves, the color-octet mechanism is not essential for perturbative consistency, 
but it may be important nonetheless in cases where the color-singlet mechanism is suppressed 
by a short-distance coefficient. A global analysis of the data on ip' production from all high 
energy processes could confirm the importance of color-octet mechanisms in the production 



of the S-wave states of heavy quarkonium. 
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Figure Captions 

1. CDF data for prompt if)' production compared with theoretical calculations of the 
total leading order contributions (dotted curves), the total color-singlet fragmentation 
contributions (dashed curves), and the color-octet fragmentation contribution (solid 
curves). The normalization of the color-octet fragmentation contribution has been 
adjusted to fit the data. 
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